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1  INTRODUCTION 


Inclusion  of  the  effects  of  clouds  in  climate  models  to  a  remarkably  high 
degree  of  precision  is  almost  certain  to  be  a  prerequisite  to  any  reasonable 
degree  of  reliability  in  climate  prediction  because  of  the  sensitivity  of,  for 
example,  surface  temperatures  to  cloud  albedo.  Unfortunately,  many  of  the 
majoi  physical  processes  taking  place  in  clouds  are  still  poorly  understood. 
Some  of  them  involve  not  only  physics,  but  marine  biology,  oceanic  and 
atmospheric  chemistry,  and  other  disciplines  outside  of  physics,  as  well.  It  is 
likely,  therefore,  that  not  only  do  physical  processes  that  are  understood  need 
to  be  modeled  much  more  accurately  than  is  now  done,  but  many  presently 
ignored  processes  have  to  be  included  as  well.  Indeed,  this  is  the  motivation 
behind  DOE’s  Atmospheric  Radiation  Measurement  (ARM)  program. 

This  paper  addresses,  as  an  illustrative  example,  one  such  previously 
unknown  complex  interdisciplinary  process  providing  a  feedback  loop  which 
may  have  a  major  impact  on  the  effect  on  global  climate  of  the  steadily 
increasing  growth  of  greenhouse  gases  in  the  atmosphere.  It  is  to  be  stressed 
that  this  is  only  one  of  a  number  of  such  feedback  loops,  many  of  which  have 
probably  not  even  been  thought  of  yet,  but  all  of  which  are  entirely  ignored 
in  present  day  computer  models. 

Clouds  influence  climate  in  two  major  ways;  they  reflect  incident  short 
wavelength  solar  radiation  thus  preventing  it  from  reaching  the  earth’s  sur¬ 
face  and  ;;  atiag  it,  and  they  absorb  outgoing  long  wavelength  radiation  from 
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the  earth  thus  reducing  the  earth’s  ability  to  cool  itself.  Which  of  these  two 
competing,  opposite  sign,  effects  dominates  is  a  sensitive  function  of  the  in¬ 
teraction  of  clouds  with  radiation,  which  is  itself  a  sensitive  function  of  the 
processes  going  on  within  various  kinds  of  clouds. 

Solar  radiation  is  (to  a  good  approximation)  a  black-body  spectrum  at 
5,780°K,  whitii  peaks  in  the  visible  at  a  wavelength  of  about  0.5  fim.  At  the 
earth,  the  flux  of  solar  radiation  is  1,370  watts/m^,  and  due  to  the  existence 
of  day  and  night  and  because  the  earth  is  a  sphere,  the  average  solar  energy 
incident  at  the  top  of  the  atmosphere  is  one  fourth  of  this.  The  average  earth 
albedo  is  .3,  so  30%  of  the  incident  radiation  is  reflected,  leaving  240  watts 
m^  to  be  absorbed  by  the  earth’s  surface  and  atmosphere. 

Radiation  emitted  from  the  earth  is  (to  a  relatively  poor  approximation) 
a  black-body  spectrum  at  255'’K,  which  peaks  in  the  infrared  at  about  10  ^m. 
There  is  essentially  no  overlap  between  the  solar  spectrum  and  the  earth's 
spectrum.  As  the  outgoing  IR  radiation  pa.sses  through  the  atmosphere,  some 
of  it  is  absorbed  and  reradiated  back  tow’ard  the  earth,  but  of  course  the  net 
IR  radiation  escaping  into  space  from  the  top  of  the  atmosphere  is  also  240 
watts/m^,  since  the  earth  is  (essentially)  in  equilibrium. 

The  black-body  spectra  and  the  principal  atmospheric  absorption  bands 
are  shown  in  F  igure  1 . 

The  albedo  of  the  earth’s  surface  (insofar  as  it  can  be  measured  from 
satellites)  is  no  more  than  .15.  The  average  cloud  cover  over  the  earth  is 
observed  to  be  60%  (65%  over  th«»  ocean  and  50%  over  the  land).  Hence,  to 
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Figure  1.  Curves  of  black-bodY  energy  B,  at  wavelength  A  for  5,780  K  (approximating  to  the 
sun  s  temperature)  and  255  K  (approximating  to  the  atmosphere's  mean 
temperature)  The  curves  have  been  drawn  c*' equal  areas  since  integrated  over  the 
earth’s  surface  and  all  angles  the  solar  and  terrestrial  fluxes  are  equal 
Absorption  by  atmospheric  gases  for  a  clear  vertical  column  of  atmosphere.  The 
positions  of  the  absorption  bands  of  the  mam  constituents  are  marked.  (From 
R.  M.  Goudy  "Atmospheric  Radiation. "  Oxford.  1964  ) 
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make  up  the  average  albedo  of  ,3,  the  albedo  of  clouds  is  about  ,15.  Clouds 
are  therefore  the  most  important  component  in  the  amount  of  reflected  sun¬ 
light,  and  their  existence  is  crucial  in  determining  the  surface  temperature 
of  the  earth.  In  fact,  (if  other  parameters  were  held  constant)  a  change  of 
cloud  albedo  by  2%  would  warm  the  earth’s  surface  by  PC'.  Evidently,  the 
modeling  of  clouds  in  climate  models  must  be  done  very  accurately. 

Clouds  are  composed  of  water  droplets  and/or  ice  crystals,  which  form 
on  cloud  condensation  nuclei  (CCN).  (In  principle,  condensation  of  water 
vapor  can  occur  at  humidities  above  100%,  but  because  the  vapor  pressure 
increases  with  the  curvature  of  the  surface  at  which  condensation  occurs,  in 
practice  humidities  of  above  300%  are  needed  for  pure  water  vapor  to  con¬ 
dense.  Therefore  actual  cloud  formation  takes  place  because  of  the  presence 
of  relatively  large  aerosols.  These  aerosols  constitute  CC.N'.)  Clouds  occur  in 
several  different  forms,  depending  on  their  altitude,  formation  mechanisms, 
hydrodynamic  properties,  etc.  Table  shows  the  distribution  of  the  various 
forms  over  land  and  ocean.  The  relative  amounts  of  these  types  often  have 
diurnal  variations  which  are  of  considerable  importance.  For  example,  the 
ratio  of  stratus  to  stratocumulus  varies  from  less  than  ‘20%  at  local  noon  to 
nearly  30%  in  the  early  morning,  while  cumulonimbus  peaks  at  30%  in  late 
afternoon  or  evening  and  vanishes  early  in  the  morning"'. 
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Table  1 

Average  %  of  Sky  Cover  over  { 

Land 

of  Various  Cloud  Types 


2  INCIDENT  SOLAR  RADIATION 


For  incident  solar  radiation,  emission  is  irrelevant  because  at  A  ~  0.5/im. 
the  black  body  spectrum  at  the  earth’s  temperature  of  around  270°K  is  es¬ 
sentially  zero.  So  only  scattering  and  absorption  count.  To  a  first  approx¬ 
imation,  let  us  neglect  scattering  into  the  beam,  and  assume  the  incident 
solar  radiation  suffers  only  losses  due  to  clouds.  Thus  the  radiation  intensity 
I\{k,z)  of  photons  of  wavelength  A  travelling  in  the  direction  k  satisfies^ 

dh(kz)  =  ~{a^^Hz)+a^^^^Hz))^(z)Ix{k.z)dz,  (2-1) 

where  cr  are  absorption  and  scattering  cross  sections  and  .V  is  particle  density. 
Therefore, 

h{kz,)  =  Ix{k,z,)€^'^^-^'K  (2-2) 

where  the  (dimensionless)  optical  depth  t  is 

t{z2,z)^  r  crl{z)N{z)dz  (2-:l) 

•'»i 

and 

<7'  =  (2  -  4) 

is  the  extinction  cross  section. 

Each  type  of  cloud  has,  at  a  given  height  :rj,  a  distribution  n(r,  z)  of  sizes 
with  particles  of  radius  r  (ice  crystals  are  of  course  not  spherically  symmetric, 
though  their  absorption  cross  section  does  not  differ  markedly  from  that  of 
liquid  water  droplets).  In  general,  the  cross  section  (t  will  also  depend  on  the 
particle  size.  For  wavelengths  that  are  small  compared  to  the  particle  size, 
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as  is  typically  the  case  for  solar  radiation  (see  Figure  2),  the  cross  section 
is  twice  geometrical:  cr  =  2rrr^.  Thus  (if  we  take  2,  to  be  the  top  of  the 
atmosphere  where  7  =  0)  the  optical  depth  at  height  r  is 

/•top  r 

7(5)  =  27r  /  dz'  I  drr^n(r,z').  (2  —  5) 

and  n(r,z')  is  the  droplet  size  distribution,  so  that  N{z')  =  zn{r,  z')dz' . 
This  equation  is  usually  re-expressed  in  terms  of  the  liquid  water  content 
(LWC)  of  the  cloud,  defined  by 

L\VC{z)  =  p,„  J  n(r,z)dr.  12  -  6) 


where  ■=  lO^’p/rrr^  is  the  density  of  water.  The  effective  radius  is  define 
by 

LWC(z)  -It 


Pu- 


)dz 


'eff'-') 

T  herefore  one  finally  writes  the  optical  depth  in  the  form 

T  I 

pu 

If  is  below  the  cloud,  then  the  total  optica!  thickness  of  the  cloud  is 

.  3  lAVP 


(2  -  7) 


■(')  =  ^  /  L\VC(c')d7'. 

2  Pu  r„(r  J- 


(2  -  8) 


;'ff 


(2  -  9) 


where  the  liquid  water  path  is 


rtop 


lAVP  =  /  ■  lAVC(c')  dz' 
y  bottom 


(2-10) 


(’rudely,  then.  i\V  P  -  LWC  ■  /,  where  i  is  the  cloud  thickness  The  liquid 
water  content  varies  greatly  with  cloud  type,  as  shown  in  Table  2'*,  as  does 
the  mean  cloud  thicknes.s.  Finally,  when  al!  of  this  is  put  together,  and 
the  average  solar  zenith  angle  of  GO®  is  includf'd,  the  previously  inferred 
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Table  2 


LWC  (g/m^) 

t  (km) 

Cumulus 

.4 

2 

Cumulonimbus 

Tropical 

1.0 

5 

Trade  Wind 

1.5 

2 

Midlatitude 

1.5 

2.5 

Polar 

1.5 

2 

■•Mtostratus/ Altocumulus 

.1 

.5 

Strat  u.s/Stratocumulus 

.2 

1 

Nimbostratus 

.1 

3 

average  cloud  albedo  of  .45  results  from  an  average  optical  thickness  t*  =  6, 
a.ssuming  a  simple  geometry  of  plane  parallel  clouds.  Since  in  fact  clouds  are 
horizontally  variable,  and  there  is  some  absorption,  this  value  is  actually  a 
lower  bound. 

Experimentally  measured  cloud  albedos  vary  greatly,  as  a  function  of 
type,  latitude,  solar  zenith  angle,  and  other  parameters.  But  overall,  they 
are  not  inconsistent  with  the  inferred  value  of  .45.  Therefore,  we  may  have 
some  confidence  in  the  value  of  r*  obtained  above. 
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From  tlu'  formula  for  r*,  we  set'  that  it  varies  inversely  with  liie  eire<’tiv(‘ 
'Iroplet  radius,  rims  a  cloud  having  the  same  lAVP  as  another  but  a  large 
number  of  smaller  droplets,  will  have  a  larger  r*  and  tlierefore  a  larger  albedo. 

If  one  asks,  therefore,  what  changes  may  take  [)lace  in  cloud  albedo 
(to  which,  as  we've  remarked  earlier,  the  earth's  surface  temperature  is  <'x 
tremely  sensitive)  due  to  man's  activiti(*s,  we  must  concentrate  on  r,ji-  For 
example,  an  increase  in  j\jf  of  10%  (which  corresponds  to  a  decrease  in  A.N 
of  30%)  reduces  the  albedo  enough  so  that  A7'  surface  =  1.3°('.  Or.  as  an 
other  example,  reducing  ,V  by  a  factor  of  2  is  ot|uivalent  to  <ioubling  the  CO.; 
concentration  in  the  atmosphere. 

A  further  effect  of  reducing  droplet  size,  while  krvping  the  total  L\\  ( ' 
constant,  is  likely  to  be  an  iroTease  in  cloud  lifetime,  and  conse'quently  of 
average  cloud  cover.  Figure  2’  shows  that  typical  droplet  sizes  are  now  .>10 
/mi.  't  hese  do  not  rain  out  until  they  coalesce  to  form  larger  dro|)s  of  -oO  //m 
diameter  or  more.  I  herefore  if  the  mean  droph't  size  were  to  decrease,  the 
tim<'  to  coalesce  to  a  size  which  rain  will  grow,  thus  iiuri'asing  average  cloufl 
lifet  ime. 

,V  and  r,.jj  are  largely  determined  by  the  number  of  (’(’N  available. 
These  vary  widely  between  land  and  water  (Figiiri*  3)'’. 

Overall,  from  the  foregoing  discussion  of  <  loud  aliiedo.  any  model  pur¬ 
porting  to  predict  global  climate  cliaiige  <lue  to  greenhouse  gas  forcing,  or 
anything  else,  rn-ecJs  to  be  able  to  evaluate  the  fractional  surface  areas  of  the 
earth  covered  by  (particularly  low)  douds,  the  liquid  water  content  of  the 
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Figure  2.  Average  drop  size  spectrum  for  the  arctic  stratus  clouds 
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Figure  3.  (a)  Percentage  of  marine  cumulus  clouds  with  indicated  droplet  concentrations  (b) 

Droplet  size  distributions  iti  a  marine  cumulus  cloud  (c)  Percentage  of  continental 
cumulus  clouds  with  indicated  droplet  concentrations  (d)  Droplet  size  distributions 
in  a  continental  cumulus  cloud  Note  change  in  ordinate  from  (b) 
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clouds,  the  droplet  effective  radius,  the  droplet  number,  and  the  number  of 
CCN,  to  an  extraordinarily  high  degree  of  precision;  less  than  5%  accuracy 
will  be  required  in  all  of  these  quantities. 

This  is  an  exceedingly  stringent  requirement  on  GCMs,  and  on  comput¬ 
ing  capacity.  Nothing  like  this  precision  is  now  available,  and,  indeed,  many 
of  these  parameters  are  not  even  included  in  present  models. 
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3  CLOUD  CONDENSATION  NUCLEI  OVER 
THE  OCEAN 


As  greenhouse  gases  in  the  atmosphere  increase,  the  earth  warms,  evap¬ 
oration  increases,  and  the  liquid  water  content  of  all  clouds  except  cirrus 
increases.  Therefore,  both  the  albedo  and  the  absorption  of  outgoing  IR 
radiation  increases.  Most  GCMs  predict  that  the  net  effect  is  a  positive  feed¬ 
back,  because  of  the  simple  observation  that  an  increase  of  temperature  at 
low  altitude  will  decrease  low  clouds  while  a  decrease  of  temperature  at  high 
altitude  will  increase  high  clouds.  Various  effects,  however,  may  reverse  this. 

The  impact  of  cirrus  clouds  is  one  uncertainty.  Their  effect  is  diffi¬ 
cult  to  compute,  since  they  are  composed  entirely  of  ice  crystals,  which  are 
anisotropic  and  whose  effect  on  radiation  is  not  quantitatively  well  under¬ 
stood.  They  are  also  not  well  studied  experimentally;  they  are  very  high  and 
also  often  even  difficult  to  see  visually. 

It  is  also  unlikely  that  fractional  cloud  cover  will  remain  unchanged  if 
evaporation  increases.  Generally,  an  increase  in  cloud  cover  will  be  a  negative 
feedback.  Figure  4^  shows  measured  changes  in  average  cloud  cover  over  the 
ocean  from  1952  and  1980,  as  a  function  of  latitude,  annually  averaged.  Could 
this  be  due  to  global  warming? 

A  major  uncertainty  is  the  effect  of  global  warming  on  the  number  of 
CCN,  particularly  oceanic  CCN. 
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Figure  4.  Change  in  %  •  cloud-cover  from  1952  to  1981  (Annual) 
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Marine  CCN  have  two  major  components;  sea  salt  and  non-sea  salt 
sulfate  (abbreviated  NSS).  By  volume,  these  are  comparable.  But  the  sea 
salt  component  is  composed  of  much  larger  particles,  so  by  number  NSS 
completely  dominates  (see  Figure  5)®. 

To  act  as  a  CCN,  an  aerosol  particle  must  be  hydrophilic  and  above 
a  critical  size,  which  is  a  function  of  the  degree  of  supersaturation  of  water 
vapor  in  the  cloud.  Over  the  ocean  the  critical  size  is  thought  to  be  in  the 
range  of  .05  to  .14  pm,  corresponding  to  supensaturation  of  0.1%  to  0.5%. 
(See  Figure  5  again.) 

Empirically,  the  number  of  CCN  over  the  ocean  is  around  100  per  cm 
(to  be  compared  with  tens  of  thousands  per  cm^  over  polluted  land  areas), 
and  since  this  is  about  the  same  as  the  droplet  number  density  in  marine 
clouds,  it  is  thought  that  the  number  of  CCN  available  is  a  limiting  factor  in 
the  growth  of  marine  clouds. 

If  we  accept  the  idea  that  NSS  are  the  dominant  CCN,  we  must  next  ask 
how  a  change  in  global  climate  will  affect  the  number  of  NSS  particles.  (We 
recall  the  apparently  very  strong  correlation  with  the  ice  age  ending  15,000 
years  ago  shown  in  Figure  6  )® 

The  present  concept  is  that  NSS  is  produced  by  the  oxidation  of  various 
sulfur  gases  coming  from  the  surface  of  the  sea®.  Dimethyl  sulfide  ((CH3)2S, 
abbreviated  DMS)  is  alleged  to  be  the  major  oceanic  source,  and  it  is  be¬ 
lieved  to  be  of  biologic  origin.  It  has  a  measured  concentration  of  about 
lOOng/liter  at  the  sea  surface,  which  varies  relatively  little  over  the  whole 
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(a)  Marine  Aerosol  Volume  Distribution 


(b)  Marine  Aerosol  Number  Distrftxjtion 


Figure  5. 
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ocean  surface  (by  a  factor  of  no  more  than  about  2).  The  entire  chain  of 
marine  biochemistry  leading  to  this  concentration  of  DMS  is  very  complex, 
involving  a  number  of  other  molecules  and  many  types  of  marine  animals 
and  plants.  In  Figure  7®,  displayed  is  the  presently  conjectured  life  cycle  of 
D.MS  in  the  sea. 

For  our  purposes  it  is  not  necessary  to  understand  this  intricate  network 
of  connections.  The  critical  question  is  whether  or  not  oceanic  warming  would 
increase  or  decrease  DMS  production,  and  with  it,  the  number  of  CC.N  and 
hence  both  the  albedo  and  cloud  lifetime. 

Figure  6  suggests  a  positive  feedback;  a  colder  environment  produces  an 
increase  in  sulfate  CCN,  which  in  turn  increased  cloud  albedo  which  further 
cooled  the  earth.  But  in  fact  which  way  DMS  production  will  change  if  the 
earth  were  to  warm  is  not  yet  known.  The  effect  could  be  either  a  positive  or 
a  negative  feedback,  and,  given  the  extreme  sensitivity  of  the  earth’s  surface 
temperature  to  cloud  albedo,  as  discussed  earlier,  the  feedback  effect  could 
be  very  important. 

Evidently,  further  research  is  needed  on  this  topic,  and  very  likely  it 
must  be  incorporated  in  some  way  or  other  in  computer  models. 

Finally,  the  existence  of  this  potentially  large  feedback  mechanism,  so 
recently  discovered,  suggests  that  there  are  other  important  feedbacks  which 
are  still  unnoticed.  We  are  not  yet  ready  for  computers  to  hand  us  believable 
predictions  about  future  climates,  even  on  very  gross  scales. 
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Figure  7. 
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